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Abstract 

We study the generalized parton distributions, including the helicity-flip ones, 
using Nambu-Jona-Lasinio model based on a relativistic Faddeev approach with 
'static approximation'. Sum rules relating the generalized parton distributions to 
nucleon electromagnetic form factors are satisfied. Moreover, quark-antiquark con- 
tributions in the region — £ < x < £ are non- vanishing. Our results are qualitatively 
similar to those calculated with Radyushkin's double distribution ansatz using for- 
ward parton distribution functions calculated in the NJL model as inputs. 
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1 Introduction 



Together with static properties, electromagnetic form factors and parton distribution 
functions are traditionally the main sources of information on the internal structure of 
the nucleon. The electromagnetic form factors of the nucleon describe the distributions of 
charge and magnetization within the nucleon, and they are determined from the electron- 
nucleon elastic scattering. Nucleon parton distribution functions are measured in deep 
inelastic scattering of leptons. It is well known that in the Bjorken limit, the deep in- 
elastic scattering data can be interpreted in a simple and intuitive picture of incident 
letpons scattered by point-like and asymptotically free partons inside the nucleons. The 
parton densities extracted from these processes encode the distributions of longitudinal 
momentum and polarization carried by quarks, antiquarks, and gluons within a fast mov- 
ing nucleon. 

With the advent of a new generation of high energy high luminosity lepton accelerators, 
a wide variety of exclusive processes in the Bjorken limit become experimentally feasible 
[1, 2]. Theoretically, it has been shown that [3, 4, 5], just like deep inelastic scattering, 
these exclusive processes are also factorizable within the framework of perturbative QCD, 
so that the hard (short-distance) part is calculable, and the soft (long-distance) part 
can be parameterized as universal generalized parton distributions (GPDs). The GPDs 
provide information on parton transverse as well as longitudinal momentum distributions. 
Furthermore, apart from the parton helicities, they also contain information on their 
orbital angular momenta. Hence measurement of GPDs would allow us to determine the 
quark orbital angular momentum contribution to the proton spin, and provide a test of the 
angular momentum sum rule for proton [6] . In the limit of vanishing transverse momentum 
transfer, the GPDs reduce to the familiar parton distribution functions. Furthermore 
their first moments give the nucleon form factors [7]. So the GPDs provide a connection 
between nucleon properties obtained in inclusive (parton distributions) and elastic (form 
factors) reactions, giving us considerable amount of new information on the structure of 
the nucleon. Excellent reviews on the GPDs can be found in Refs. [8, 9, 10]. 

GPDs can be measured in deeply virtual Compton scattering (DVCS) and in hard 
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exclusive leptoproduction of mesons. First measurements of DVCS related to GPSs have 
been recently reported in [13, 14, 15], and other experiments designed to measure GPDs 
in exclusive reactions are expected to be carried in the near future [1, 2, 16, 17]. Hence 
theoretical estimates of the GPDs will provide a very useful guide to future experimental 
efforts. 

Unlike parton distributions, the GPDs in general cannot be interpreted as particle 
densities, they are instead probability amplitudes. However, like the parton distribution 
functions, GPDs reflect the low energy internal structure of the nucleon, and are at present 
not directly calculable from first principle in QCD. A first attempt to calculate the first 
moments of GPDs in quenched lattice QCD at large pion mass has recently been reported 
[18]. However, there is still a considerable large gap in quark mass to bridge between the 
state-of-art lattice QCD calculations and the chiral limit. Other theoretical calculations 
have also been performed in various QCD-motivated models of hadron structures such as 
MIT bag models [19, 20], chiral quark-soliton model [21], light-front model [22], Bethe- 
Salpeter approach [23], and constituent quark models [25, 26]. In this work, we calculate 
the nucleon GPDs in the NJL model [28]. 

One of the most important features of QCD is chiral symmetry and its spontaneous 
breaking which dictate the hadronic physics at low energy. As an effective quark theory 
in low energy region, NJL model [28] is known to conveniently incorporate these essential 
aspects of QCD. Models based on the NJL type of Lagrangians have been very successful 
in describing low-energy mesonic physics [29] . Based on relativistic Faddeev equation the 
NJL model has also been applied to the baryon systems [30, 31]. It has been shown that, 
using the quark-diquark approximation, one can explain the nucleon static properties rea- 
sonably well [32, 33]. If one further takes the static quark exchange kernel approximation, 
the Faddeev equation can be solved analytically. The resulting foward parton distribu- 
tion functions [34] successfully reproduce the qualitative features of the empirical valence 
quark distribution [35] . Recently, NJL model has been used to investigate the quark light 
cone momentum distributions in nuclear matter and the structure function of a bound 
nucleon [36] as well. In this work, we extend such a NJL-Faddeev approach to calculate 
the nucleon GPDs. Since NJL model is a relativistic field theory, the GPDs so obtained 
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will automatically satisfy all the general properties such as the positivity constraints and 
sum rules [10]. 

This paper is organized as follows: in Section 2 we explain the model used in this 
work. In Section 3 we outline the calculation of GPDs. Results and discussions are given 
in Section 4, and finally a summary is given in Section 5. 

2 The NJL model for the nucleon 

The SU(2)j NJL model is characterized by a chirally symmetric four-fermi contact inter- 
action Lagrangian Cj. With the use of Fierz transformations, the original NJL interaction 
Lagrangian Ci can be rewritten in a form where the interaction strength in any channel 
can be read off directly [31]. In particular, we are interested in the following channels: 

Ci,« = G n [{^f - (^ 75 r^) 2 ] , (2.1) 
Ci,s = G s [4,{^C)T 2 p A i, T ] [^(C-^r^} , (2.2) 

where f3 A = \J\\ A (A=2,5,7) are the color 3 matrices, and C = 27270- £i,n represents 
the interaction in the + and 0~ qq channels corresponding to the sigma meson and the 
pion, respectively. Cj jS describes the qq interaction in the scalar diquark channel (J 71- = 
+ ,T = 0). The interactions (2.1) and (2.2) are invariant under chiral SU(2)l x SU(2)r 
transformation. The coupling constants G n and G s are related to the ones appearing in 
the original Ci by Fierz transformation. Here we shall for simplicity take £/ )7r and Ci tS 
as starting points, and treat and G s as free parameters. 

The reduced t-matrices in the pionic and scalar diquark channels are given by the 
following expressions [34]: 

rAk) = i + 2(?!n[(F)' TD(k) = i + 2Gji D (kiy (2 ' 3) 

with the "bubble graph" contribution given by 

Tl^k 2 ) = Tl D (k 2 ) =GiJ 0- itlD [ l5 S(q) l5 S(k + q)], (2.4) 

where S(q) = 1/0/— Mq + ie) is the Feynmann propagator and Mq is the constituent 
quark mass. 



In order to simplify the numerical calculations we approximate Tn(k) by 

ro(k) - «G, - (2.5) 

in the actual calculation, where M D is the diquark bound state mass. g|, is the residue 
of the pole of Tn(k), 



fc 2 =M£, 



In performing the four-momentum loop integral, we have to introduce a cutoff scheme. 
In this paper we will adopt the Pauli-Villars (PV) regularization scheme [37] which 
preserves the gauge invariance and can be applied to light-cone (LC), Euclidean, and 
Minkowsky space integrals. We will follow [38] to determine the subtracted terms in PV 
regularization scheme. 

The original definition of PV regularization scheme is defined by the following substi- 
tution in every loop integral: 

1 1 " [ 1 1 \ 

k\ - M 2 ' ' ' k% - M 2 f^ Ci \ k\ - M 2 - A? ' ' ' k% - M 2 - A? J ' ( } 

where c = 1 and A = 0. For the convergence of the loop integrals, we need to impose 
at least 2 conditions 

£> = 0, X>A 2 = 0. (2.8) 

i=0 i=0 

Thus, we need at least 2 subtractions. In order to reduce the number of parameters, 
we choose n = 2 and take the limit A 1 — > A 2 = A. The reduction formula for PV 
regularization scheme then becomes 

E cJ(A') = /(0) - /(A 2 ) + A 2 ^^. (2.9) 
i=o aA 

In Ref . [34] , the relativistic Faddeev equation was solved analytically under the static 
approximation, namely the momentum dependence of the quark exchange kernel is ne- 
glected, i.e., 

1-1 . . 

(2.10) 



tf g -M Q + ie M Q 
The analytical solution for the quark-diquark T-matrix is given by 



T( p ) = — i — , (2.11) 
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where Hn(p) is the quark- diquark bubble: 

n*(p) = f-0^S(k)r D (p-k). (2.12) 

The nucleon mass M N is obtained from the pole of quark-diquark T-matrix of Eq. 
(2.11), whose behavior near the pole is given by 

T{p) - £ T n(p, s)T N (p, s)/(p 2 -M 2 n + ie), (2.13) 

s 

where f N = rjv7 . 

Together with Eq. (2.11), it leads to the following expression for the nucleon vertex 
function T N (p, s): 

T N (p,s) = y/z^u N (p,s), (2.14) 
dU N (p)Y 1 



Z 



dp 1 



(2.15) 

i>=u N 



where un(p,s) is the nucleon Dirac spinor with normalization un(p, s)un(p, s) = 2M^. 
With this normalization convention, un(p, s)'J ± un(p, s) = 2p ± and the nucleon vertex 
satisfies the relation 

_1 -f N {p, s)^^-T N ip, s) = 1, (2.16) 



2p_ dp + 

where the LC variables are defined by a ± = a T = (a ± a 3 )/V2, and a± :i = — a? ± for 
i = 1,2. 

3 GPDs of the nucleon 

It is well known that inclusive deep inelastic scattering of leptons from nucleon is described 
by universal parton distribution functions. Hard exclusive processes measure another kind 
of structure functions called generalized parton distributions (GPDs) of the nucleon; they 
can be diagrammatically depicted in Fig. 1. Just like ordinary parton distributions, the 
GPDs are also process independent universal functions. 

In hard exclusive processes, a high energy virtual photon of momentum g M is absorbed 
by a quark in a nucleon, producing a real photon or a meson and without breaking up 
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k=K-A/2 y k 



k'=K+A/2 



p=P-A/2 



\ 



p'=P+A/2 

Figure 1: Soft amplitude for the GPDs. 



the nucleon [8, 9, 10]. It is customary to choose a frame where the averaged nucleon 
four-momentum P = (p + p')/2 and are collinear along the z-axis [4]. Then the GPDs 
H(x, £, A 2 ) and E(x, £, A 2 ) are formally given by the leading twist (twist-two) part of the 
following amplitude 

2k 



J dy + e* xP -y+ < p'A'|^(-y/2)7 + ^%/2)bA > 



y-=y±=o 



u N (p', A') 



2M, 



N 



u N (p,X) + 



where A = p' —p, superscript q denotes the quark flavor, \pX > stands for a nucleon state 
with momentum p and helicity A, and the meaning of x and £ will be made clear in the 
momentum representation below. The ellpsis (• • •) denotes the higher-twist contributions. 
In momentum space the above expression can be written as 



u>n(p', A') 



2M, 



N 



un(p, A), 



= J0^A x - K+ / P+ ) tr h + X q N(p, P ',K)} (3.17) 

where k — (x + £)-P + and k' — (x — £)P + are respectively the initial and final quark 
momenta, K = (k+k')/2, and XqN&p' , K)^ = J d 4 ye iK ' y < p'\'\i> i (-y/2)^ j (y/2)\p\ > is 
the quark-nucleon scattering amplitude. The LC momentum fraction x and the skewness 
£ are given by x = K + / P + and £ = -A+/(2P + ), with 



0<£< 



-A 2 



AM% - A 2 



< 1. 
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Due to the on-shell conditions, p = p 

A 2 = 



= Mjv, we also have: 
IC'.Uv + A 2 



(3.18) 



(3.19) 



It is then clear that a GPD describes the amplitude of emitting a parton with mo- 
mentum fraction x + £ in a nucleon and reabsorbing one with momentum fraction x — £. 
If x > £, both the emitted and absorbed partons are quarks; if x < — £ then both are 
antiquaries. Finally, if \x\ < £, the two partons involved are a quark- qntiquark pair. From 
this physical picture, it is clear that, in the forward scattering limit £ = 0, GPDs reduce 
back to the familiar parton distributions q(x): 



q(x) = H q (x,0,0). 



(3.20) 



Furthermore, by integrating H q (x, £, A 2 ) over x, we recover the nucleon elastic form fac- 
tors: 

-l 



J i dxH0(x,£,A 2 ) = F?(A 2 ) 
| 1 i ^(x,£,A 2 ) = F|(A 2 ). 



(3.21) 
(3.22) 



In the NJL model, the GPDs can be calculated by evaluating the Feynman dia- 
grams shown in Fig. 2, where the contributions from the quark and diquark currents, 
Jy A (x, £, A 2 ) and Jy \(x, £, A 2 ), are shown separately. Note that in the NJL model we 
use here, only the isoscalar diquark is considered, then it is easy to see that: 



J£>,£,A 2 ) = J« A (x,£,A 2 ) + Ji? iA (a;,£,A 2 ), 
4,a(^,A 2 ) = Ji?, A (x,£,A 2 ), 



(3.23) 
(3.24) 



where superscript Q(D) denotes the quark (diquark) current contribution. We further 
write (X = Q, D) 



J X tX (x,C,A 2 ) = u N {p',X') 



H x (x,t^h + + E x (x,t^) ia+ ' /X 



1M, 



N 



u N (p,\). (3.25) 



Using the table of matrix elements listed in Appendix A, we can separate the left hand 
side of Eq. (3.17) into helicity conserving and helicity flipping contributions: 



M, 



N 



-u N (p', X')u N (p, A) 



x 



5y,x ((1 - e)H x (x, £, A 2 ) - eE x (x, £, A 2 )) - 5y^ x E x (x, £, A 2 ) 



(3.26) 




Figure 2: Graphical representation of the quark GPDs of the nucleon. The single (double) 
line denotes the constituent quark propagator (diquark t-matrix). The operator insertion 
stands for r y + S(K^ — xP){\ ±r z )/2 for the u(d) quark. Initial (final) nucleon momentum 
and helicity are denoted as p (p') and A (A'), and the four- momentum transfer is given by 

In the following we shall only give an outline of the calculations, and leave the details 
to Appendices B and C. 

Using simple Feynmann rules, we can directly read off the quark current contribution 
J« A (x,e,A 2 ) from Fig. 2, 

J^ x (x,Z,A 2 ) = -Z N u N (p>,\>) J 0^s{x-^j S(k')^S(k)r D (p-k)u N (p,\), 

(3.27) 

where td is the reduced t-matrix of the diquark. td can be decomposed into two terms: 

t d = t% + t£ (3.28) 

where r£ and are respectively the "contact" and "pole" contributions, as given in 
Eq. (2.5) Accordingly, the quark current contribution J Q can also be separated into two 
terms: 

J?, tX (x,Z, A 2 ) = 0(-C <x< £)J$'£(x, £, A 2 ) + 0(-£ <x< i)J$'£(x,£,A 2 ), (3.29) 

where we see that J® ,c contributes only in the region — £ < x < £, while J®' p only in 
— £ < x < 1. Thus we see that, unlike other calculations using non-relativistic quark 
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models, the field theoretic NJL model gives a non-zero contribution in the — £ < x < 1 
region. 

Our final expressions for the Hq and Eq are given in Eqs. (B.56) and (B.57). 
Similarly, we can also write down the diquark current contribution, 

Ji?, A (x,£,A 2 ) = -Z N u(p>,\>) J ^-iS(P-T)r D (t')r D (t) 



x i 



-tr 



^Cr^Sik'h+Si^C-^^ASit - k) T 5 [x - u(p, A), 



K_ 



(3.30) 



where t = T — A/2,t' = T + A/2 are the diquark momenta. 
We define two additional LC momentum fractions y, z, 



so that 



y = K + /T + , z = T + /P + , 



dy I dz5 (y - |-j s(z- ^ j = 1. 



(3.31) 



(3.32) 



Inserting Eq. (3.32) into Eq. (3.30), we can rewrite the diquark current contribution in a 
convolution form: 

Jvtfrt,* 2 ) =j0^J d yj dz5(x-yz)^/ x N (z,£,T,A) 
x (F s D (y,(,T,A)T + + F a D (y,(,T,A)A + ), (3.33) 



with 



(F s D (y, C, T, A)T+ + F®(y, (, T, A) A" 



= Wd 



d A K 



S[y 



T+ 



tr 



^CT 2 p A S^) 1 + S{k)C- 1 1 & T 2 p A S{t - kf 



(3.34) 



^ (z,t,T,A) = -tg D 2 Z N u(p',X')5 (z-^j S(P - T)r D (t')r D (t)u(p, A), (3.35) 

where go is given in Eq. (2.6), and ( is the skewness defined by the relation A + = —2(T + . 
From the Ward identity for the diquark-diquark-photon vertex T^ D {T, A), 



A M T D (Or£(T, A)r D (t) = +2i[r D (t') - r D (t)], 



(3.36) 
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we obtain 

2t[r D (t')-r D (t)} 
r D (t)T D (t)= A#|r£(TjA) ■ (3-37) 

In order to reduce the complexity of the calculation, we introduce the 'on-shell diquark 
approximation', i.e., t 2 = t' 2 = Mp. Then the vertex can be expressed in terms of a 
single form factor Gf 

F%(T,A) ~ Gf (A 2 )T M . (3.38) 

(see Appendix C for details). 

Using Eqs. (2.5) and (3.38), we can rewrite Eq. (3.37) as 

T ^>^ ~+ G?(Ai)(t> 2 -M 2 D )(t 2 -M 2 D y (3 ' 39) 
Substituting this result into Eq. (3.35), we finally arrive at 

^ (z,Z,T, A) = iZ N u( P >, AO g *<* " 11^2?)% I T Ml) u ^ A ) • ( 3 - 4 °) 

The final results for H D (x, £, A 2 ) and E D (x, £, A 2 ) are given in Eq. (C.63). In Appendix 
C, apart from A M , we have introduced another momentum transfer variable A^ inside 
the convolution integral. In a complete evaluation, we should have A M = A^>. However 
the on-shell diquark approximation gives raise to an ambiguity. In [26] , it is assumed that 
A 2 = A 2 D and A + = A J, then ( — £/z. This form is adopted for small £ and A 2 ± « M%. 
However this choice of A£> is not satisfactory in our case, since for A + = Aj implies that 
the GPDs are non-vanishing only in the region £ < z, and it follows from Eq. (C.69) that 
the sum rule is explicitly broken. In order to preserve the sum rule relation, which we 
believe is important, we let A^ = A 2 DA _. Then ( is fixed by the relation A 2 = A 2 D because 

A2 \e\ii ■ a 2 



i-e (3 ' 41) 

A 2 =-!^±#k. (3.42) 



Finally we also include the photon-quark vertex correction for the GPDs and electro- 
magnetic form factors which arises from the structure of the constituent quark. Specifi- 
cally, we sum the series of ring diagrams as shown in Fig. 3. In the spirit of vector domi- 
nance, we demand that the resultant photon vertex possesses a pole at A 2 = M 2 (A 2 = 
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Mp) in the isoscalar (isovector) channel . This effectively replaces the bare vertex r a 7 M 
by 

(3.43) 



l + 2G a A2n y , T (A 2 )' 

where a = 0(i) corresponds to the isoscalar (isovector) part (r° = 1), and the correspond- 
ing coupling constants are G w (for a = 0) and G p (for a — i). The definition of Hv,t(A 2 ) 
and the details of the calculation can be found in appendix D. 



+ 




+ 




p,C0 



+ 



Figure 3: The vector meson dominance corrections to the •yqq vertex. The dotted line 
represents the vector mesons u and p. 



As we shall see in the next section, these vertex corrections significantly improve the 
momentum dependence of the electromagnetic form factors calculated in the NJL model. 

4 Results and discussions 

In this section, we present the results of our calculation. We shall first explain the choice 
of parameters in our model. Subsequently, numerical results are presented and compared 
with those obtained from other works. 

In the NJL model we have adopted here, the constituent quark mass is taken to be 
Mq = 400 MeV, which is within the range of values used in other works [40]. Using this 
constituent quark mass, together with the pion mass = 140 MeV and the pion decay 
constant f w = 93 MeV , we can determine the Pauli-Villars cutoff parameter A = 739 
MeV, the coupling constant G n = 10.42 GeV -2 , and the current quark mass m q = 9 MeV. 
Furthermore, we set G s = 0.65G7,- so that the solution of the Faddeev equation reproduces 
the experimental nucleon mass M N = 940 MeV, the scalar diquark mass is then fixed to 
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be M D = 590 MeV. Finally, the coupling constants G w = 7.34 GeV~ 2 and G p = 8.38 
GeV -2 are determined from the poles of Eq. (D.73), so that the physical vector meson 
masses m w = 783 MeV and m p = 770 MeV are reproduced. 

The nucleon electromagnetic form factors Ff (A 2 ) and F™(A 2 ) calculated in the NJL 
model are depicted in Figs. 4(a) and 4(b), respectively, where results with and without 
corrections to the photon-quark electromagnetic vertex (see Fig. 3) are shown separately. 
For comparison, we have also plotted the results corresponding to the familiar " dipole fit" 
to the experimental data of Gg's and Gm's by dashed lines: 

<? B ,G$ oc (1-A 2 /A 2 )" 2 

G n E = 0, (4.44) 

with Aq = 0.71 GeV 2 . We see that the effect of the vertex corrections is rather sizable. 
In the proton case, where the data are much more precise, inclusion of the vertex cor- 
rection significantly improves the agreement with experimental data, which is very well 
parameterized by the "dipole fit". Nevertheless discrepancies still exist for —A 2 > 0.5 
GeV 2 . Note that for A ± = ± , -A 2 = 0.5 GeV 2 corresponds to f ~ 0.35, and in this 
work we are only concerned with small £ (< 0.3). In the case of the neutron, the effect of 
the vertex correction is small. Compared with the "dipole fit" , the NJL model result for 
F™(A 2 ) is similar in magnitude, but opposite in sign. Unfortunately the available data in 
this case are scattered with large error bars, so that it is not possible to determine which 
curve fits better. 

Similarly, the nucleon form factors F<f' n (A 2 ) are plotted in Figs. 5(a) and (b). We 
see that the NJL model results are significantly different from the dipole fits in the low 
momentum region —A 2 < 1 GeV 2 . As a result, the calculated nucleon magnetic moments 
(in units of nuclear magneton) 

< JL = 1.75, ^ JL = -0.82, (4.45) 

are much smaller than the experimental values 

= 2.79, fM e n xp = -1.91. (4.46) 

This will affect the reliability of the E q (x,^, A 2 ) calculated in this model (see discussions 
below). It is known that further inclusion of the axial vector diquark channel and the 
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pion cloud are important to improve the results for the magnetic moment [34]. However, 
these effects are outside the scope of our present investigation of GPDs. 




-02 - 



-OA - 



0.5 1 1.5 2 

-A 2 [GeV 2 ] 
(b) 



Figure 4: (a) Proton form factor Ff(A 2 ). The dotted and solid lines are calculated in 
the NJL model without and without vertex corrections, respectively, and the dashed line 
is the dipole fit to the experimental data, (b) Neutron form factor ^"(A 2 ) in the same 
notation as (a). 

Having fixed the model parameters, we now present the main results of this work. The 
calculated GPD's, H u (x, £, A 2 ),H d (x, £, A 2 ), E u (x, £, A 2 ), and E d (x, £, A 2 ), are plotted in 
Figs. (6-9), for three different values of £ = 0,0.1,0.3, with A 2 given by Eq. (3.41). For 
simplicity we have assumed A^ = 0j_. 
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(b) 



Figure 5: (a) Proton form factor ^"(A 2 ), (b) Neutron form factor F^(A 2 ). Notation 
same as in Fig. 4. 
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Figure 6: H u (x,l;,A 2 ) for £ = 0,0.1,0.3. The solid lines are NJL model results, and 
the dashed lines are obtained using the Radyushkin's ansatz for the input forward quark 
distributions calculated in NJL model. 



-1 -0.5 05 1 

X 

Figure 7: H d (x,£,A 2 ) for f = 0,0.1,0.3. Notation same as in Fig. 6. 
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1 -0.5 





0.5 



Figure 8: E u (x, £, A 2 ) for £ = 0,0.1,0.3. Solid lines show NJL results while the dashed 
lines give the NJL results multiplied with a factor of F^ EX p(A 2 )/F^ Arji (A 2 ) (see text for 
explanation). The dotted line represent 




i;=0 



Figure 9: E d (x,C,,A 2 ) for £ = 0,0.1,0.3; Notation same as in Fig. 8. 
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As mentioned in Section 3, the quark-current contribution J Q to the GPDs can be 
decomposed into two terms, J® ,p and J®' c , corresponding respectively to the "pole-term" 
and "contact term" in the diquark t-matrix, see Eq. (3.29). It has been found that the 
contact term violates PCAC by as much as 13% [36] which is related to the use of "static 
approximation" for the Faddeev vertex function. Moreover, the contact term contributes 
only in the quark- ant iquark region, — £ < x < £, producing unphysical kinks at x — ±£. 
In view of these facts which indicate that the contact terms can not be assessed reliably 
in the static approximation, we have chosen to leave out the contact term contribution in 
our results. 

Like all constituent quark models, there are no intrinsic anti-quarks in the NJL model, 
therefore H Q (x, £ = 0, A 2 ) and E q (x, £ = 0, A 2 ) (q = u, d) vanish for negative x. However, 
unlike the constituent quark models, the NJL model is field theoretic in nature and the 
Fock states with antiquarks can appear in the intermediate states. Accordingly, the 
quark- ant iquark contribution to GPDs in the region — £ < x < £ is accessible in our 
calculation. 

As mentioned before, the calculated electromagnetic form factors Ff ,n (A 2 ) do not re- 
produce the experimental data in the low momentum transfer region —A 2 < 1 GeV 2 . 
These discrepancies would affect the quality of E q (x, £, A 2 ) calculated in our model be- 
cause F^A 2 ) is related to the first moment of E q (x,l;,A 2 ) through Eq. (3.22). Conse- 
quently, we scale up our calculated E q values by a factor of p^ d{ ~ exp ^> / pu,d(N r J 'L) 
them in Figs. 8 and 9. 

It is interesting to compare our results with those obtain using Radyushkin's ansatz 
[27]. Radyushkin proposed to write the GPD in terms of a "double distribution" F q ((3, a, A 2 ) 
which is assumed to be factorized: 

F q ([3, a, A 2 ) = h(a, P)q(a)F?(A 2 )/F?(0), (4-47) 

where q(x) is the forward quark distribution (or quark distribution function) and the 
profile function h(a, (5) has the property of asymptotic meson distribution amplitudes 
given in [27]: 



18 



H g (x,^, A 2 ) is then given by the convolution expression: 

H q {x, £, A 2 ) = df3 da5{x - (3 - a£)F q {(3, a, A 2 ). (4.49) 

Using the forward quark distributions calculated in the NJL model as input, we plot 
the results obtained from the above ansatz also in Fig. 6 and Fig. 7. We see that, in 
magnitudes and in shapes, Radyushkin's ansatz gives qualitatively similar results as the 
NJL model. One visible quantitative difference is that, as £ increases, the peak position 
of H q (x, £, A 2 ) shifts towards larger x in the NJL model, while it stays almost unchanged 
in the Radyushkin's ansatz. 

In Figs. (8-9), n q = F|(0). It is seen that K d d(x) is quite similar to E d (x,l; = 0, A 2 = 
0). This is because E d receives contribution only from the diquark current in our model. 
On the other hand, our result for k u u{x) is rather different from E u (x, £ = 0, A 2 = 0), in 
contrast to the results obtained with the chiral quark-soliton model [21] and constituent 
quark models [25, 26]. 

In Figs. (10-11), we compare our results with those obtained in a calculation using the 
constituent quark model [43] which is calculated using a simple gaussian wave function. 
We note that their calculation of GPDs is exactly same as [25] except the use of a different 
wave function. First of all, we see that the signs of the GPDs calculated in the two models 
agree except E u (x, £,A 2 ), that is, E u (x,^,A 2 ) calculated in the NJL model explicitly 
shows a negative contribution for small x. Secondly we see that the shifting of the peak 
positions towards larger x with increasing £ are common in both calculations. Finally, we 
observe that due to the fact that there is no quark-antiquark contribution to the GPDs in 
the constituent quark model, the curves all terminate at x — £. In contrast, as mentioned 
earlier, the NJL model is field theoretic in nature, so that quark-antiquark contributions 
is non-zero in our calculation. As a result, the range of validity is — £ < x < 1 in our 
calculation. 

Comparing our results with those obtained from the chiral quark-soliton model [21, 40], 
we find that the behavior in the range — £ < x < £ is quite different. In the case of chiral 
quark-soliton model, strong oscillatory behavior is seen around x = ±£, whereas our 
results are rather smooth. This difference arises from the fact that in the chiral quark- 
soliton model, there is a so called "d-term" contribution [9, 41, 42] which corresponds to 

19 





-1 -0.5 05 1 

X 




Figure 10: (a) H u (x,£, A 2 ) for f = 0,0.1,0.3. The solid lines are the NJL model results 
and the dashed lines are obtained with constituent quark models [43]. (b) H d (x, £, A 2 ) 
for £ = 0,0.1,0.3. Notation same as in (a). 
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Figure 11: (a) E u (x, f , A 2 ) for f = 0, 0.1, 0.3. (b) E d (x, f , A 2 ) for f = 0, 0.1, 0.3. Notation 
same as in Fig. 10. 
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the case where the active quark and antiquark are correlated in the scalar isoscalar (or a) 
channel [24]. Such a contribution is supported by the recent preliminary HERMES data 
[44] on beam-charge asymmetry in DVCS but is not included in our model. 

5 Summary 

In this work, we have calculated the spin-averaged {H q ) and helicity-flip {E q ) GPDs of 
the proton, using the NJL model based on a relativistic Faddeev approach with "static 
approximation" . The NJL model is a field theoretic approach which has been successfully 
used in the studies of the static properties and parton distribution functions of the nucleon. 
Hence the NJL model provides a reasonable framework in which to calculate the GPDs 
or off-forward parton distribution functions. Among other things, there are two major 
advantages of adopting this model. First, due to the fact that NJL model is a relativistic 
field theoretic model Fock states with anti-quarks can exist in the intermediate states, 
hence quark-antiquark contributions to the GPDs in the region — £ < x < £ are non- 
vanishing. In addition, the model independent sum rules relating the GPDs and nucleon 
electromagnetic form factors are satisfied. 

The calculated GPDs are qualitatively similar to those calculated with the Radyushkin's 
double distribution ansatz with forward parton distribution functions calculated in the 
NJL model as inputs. Comparing our results with those obtained in constituent quark 
models [25, 26], we find that the general features are similar, except for the fact that the 
region — £ < x < £ is not accessible in the latter. 

In our present treatment of the NJL model, as well as in other quark models, config- 
urations with intrinsic antiquarks are not present. Hence it is not possible to investigate 
GPDs in the region x < — £. In our case, antiquark contribution can be studied if we in- 
clude the pion cloud surrounding the three-quark core. In addition, since NJL model is an 
effective quark theory in the low energy regime, we need to evolve our results, according to 
perturbative QCD, in order to compare them with data taken in high-energy experiments. 
Such an NLO (Revolution of the calculated GPDs, from the low-momentum scale to the 
experimental one, has been carried in Refs. [26, 45]. We will leave these improvements to 
future investigations. 



22 



Acknowledgments 

The authors wish to thank T. Spitzenberg and M. Vanderhaeghen for helpful discussions. 
This work is supported, in part, by the National Science Council of ROC under grant 
Nos. NSC93-2112-M002-004, NSC93-2112-M002-058, and the Grant in Aid for Scientific 
Research of the Japanese Ministry of Education, Sports, Science and Technology, Project 
No. 16540267. 



23 



Appendices 



A Matrix elements of Dirac spinors 



Table 1 contains the matrix elements of Dirac spinors used in our calculations 
convention of [39] is adopted here: 

u N (p,X) = = {V2p + + a ± -p ± + pM N ) X (\), 



( 1 \ 



1 

72 



x(+i) 

and a±(X) and a±_ A b± are defined by 




1 

V ) 



, x(-i) 



1 

72 



/ o \ 
1 



V -i / 



1. The 



(A.50) 



(A.51) 



o_l(A) = Xa + ia , a_ L Ab ± = a b - a b . 



(A.52) 



Table 1: Matrix elements of Dirac spinors un(p' : \')M.un(p, A) 



M 


-^^-u^, X')Mu N (p, A) 


-^j^ujvG/, A')-Mujv(p, A) 
v P P 


1 

7 + 
7~ 
7 j_ • a_L 

7"7 + 
7_l • a ± 7 + 

7~7 + 7j_ • «± 
7 + 7~7j_ • «± 
a± • l±l + b± ■ 7± 


Mat _i_ Mjv 
p'+ p+ 

2 

-p± + M 2 N + i\p' ± A p ± ) 
M£ + " <Afl -L A - 


^+ (A -P± + M^ + iXp' L A p ± ) 

-Rl + «Aa ± Apx) 
2(a_|_ • 6j_ + iAa_|_ A 


P'l (A) p±(A) 
p'+ p+ 



^(p' ± (A)-p ± (A)) 
-a±(A) - ^) 

2a ± (A) 
^(p' ± (A)-p ± (A)) 

-^a±(A) 
^ai-(A) 
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B Quark current contribution 

In Eq. (3.27), the i^ + -integral can be trivially performed. Then with the help of table 1, 
we get 



u N {p\ *W + M) 7 +(# + M)u N (p, A) 



u N (p', X)u N (p, A) 



P+ 



M, 



N 



x 2 - i z )M 2 N + {2x + i- C)MM X - 



(1-x) 



+ (1 - f)(Kl + M 2 ) +(1 + x + i + e 2 )i?± • £l + iA(l + 2£ + x)K ± A Aj 
(for A' = A) 



u N (p', -X)u N (p, A) 
(for A' = -A) 



P+ 



M 



N 



2(1 - x)M N {{x + 2£^^j)M JV + Atj 



The rest of Eq. (3.27)is given by 

cPK ±J rE(p-k)+rE(p-k) 



-Z 



N 



-dK- 



(2vr) 4 (k 12 — M 2 )(k 2 — M 2 ) 
After performing the if~-integral, we obtain 

-rg(p - k) 



I 



dK 



(k' 2 — M 2 ){k 2 — M 2 ) 
-4ggg(-g < x < Q 
5Q(a;,e,^,A ± ,M 2 ,M2)' 

-Td(p - fc ) 



= F c (a;,e,^,A ± ,M 2 ,M 2 



/ 



= F P (a;,£,K ± ,A ± ,M 2 ,M 2 ) 



x 



(A;' 2 -M 2 )(A; 2 -M 2 ) 

9l 

2A Q (x, -^K ± ,-A ± ,M 2 ,M 2 D ) 

( (i-x)6(x>o + ( x + s)e(-s< x <o 



^A Q (x,^K ± ,A ± ,M 2 ,M 2 D ) B Q (x,£,K ± ,A ± ,M 2 ,M 2 D )J ' 



(B.53) 



(B.54) 



where 



Aq{x, £, K ± , Aj_, M , Mf)) = (l- x )( x -Z)(l + Z)P 2 -( x -Z)(Ki + M 2 D ) 

— * 

-(l- x )[(K L + ^) 2 + M% 
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B Q (x,Z,K ± ,A ± ,M 2 )=2 



i{x 2 - e)P 2 - i{K\ + M 2 + ^)-xK ± -A 



. (B.55) 



Combining the above results and with the help of Eq. (3.26), we arrive at the final 
expressions for H Q and E Q which can be decomposed into the 'pole' and 'contact' term 
contributions (Eq. (3.29)): H Q = H c + H P ,E Q = E c + E p . In the following we will 
explicitly write down the results of H ,p and E ,p with PV regularization scheme (Eq. 
(2.7)): 



= E 



H C ' P (x^,A 2 ) 
d 2 K, 



Ci 



A ± (A) 



M 2 N + 2x{\ - i 2 )MM N 



:i--) 2 a 



Ai + (i-e 2 )(^i+M^) + (i+x+e+m±-Aj 



x 



F c , P (x,tK ±: A ±: M 2 : M D 2 ) 

(i-a 

E c > P (x^,A 2 ) 
= E«/g^(l-*)M„ 



(B.56) 



FcAx,^K±,A ± ,M 2 ,M D 2 ), 

(B.57) 



where M 2 = M 2 + A 2 , M 2 Di = M% + A? and the q's are given in Eq. (2.7) and Eq. (2.9). 
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C Diquark current contribution 



The diquark current contribution is given in Eq. (3.34). In order to simply the calculation, 
we shall assume the initial and final diquarks are on shell, that is, t' 2 = t 2 = M|>. Then 
T ■ Ad = 0, and T 2 = M|, — A|,/4. In the above and later discussions we will explicitly 
distinguish A M and A^, since under the on-shell diquark approximation the frame where 
we calculate diquark GPDs is not necessarily the same as the one originally chosen for 
the nucleon GPDs, so that in general A^ ^ A^ D . Thus, Eq. (3.34) becomes 



{F s D (y,CA 2 D )T+ + F a D (y,(,A 2 D )A + ) 

= 191 1 WT 4 * { V ~ S~) ^ [ S Wl + S(k)S(T - K)] , (C.58) 

where A J = -2(T + , and in the frame where f ± = ± , A 2 D is given by A 2 D = _ l^a+^L . 
Integrating Eq. (C.58) over y, we can reproduce the diquark form factors G® a) 

J^dyF^(y,C,A 2 D ) = Gf ia (A 2 D ), (C.59) 

where we see that G® a (A 2 D ) is independent of ( as required. Furthermore for on-shell 
diquarks, due to the symmery under the exchange of t and t', we explicitly find that 
G?(A2,)=0. 

After integrating over K + and K~, we obtain 



(F s D (y,C,A 2 D )-2CF a D (y,C,Al))T + 
R 2 f d2K 



0(y > C) 



(2tt) 3 

C 2 (l - y) 2 T 2 + (1 - C 2 )(K 2 ± + M 2 ) + + C(l - y)K ± • Ad± 



A D (y, C, K ± , A D± , M 2 )A D (y, -(, K ± , -A D± , M 2 ) 
~ 0{\y\<C) 

(C 2 (C-y)+y 2 (i - C))T 2 + (1 + Q{Kl + M 2 ) + fl-^c-rtff, , + y K ± . a d± 

X 



A D (y, C, K ± , A D± , M 2 )B D (y, (, K ± , A D± , M 2 ) 



(C.60) 



with 

A D (yX,K ± ,A D± ,M 2 ) 
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= (V + C)(v - 1)(1 " C)T 2 + (y + ()Ki - (y - 1)(K ± - ^) 2 + (1 + C)M 2 , 
B D (y,(,K ± ,A D± ,M?) 

where T 2 = M 2 D - Ay 4. 

Next we need to calculate the following integral 

F^ x N (z,Z,T ± ,A ± ) 

i9 2 pZ N , , f^^ Siz-T+fP+WP-T) 
= ^GfW) UipA) J (t' 2 - M 2 D )(t 2 - M 2 D ) U ^ X) - 

(CM) 



Insert the PV-regularization scheme, and following the same steps as indicated in appendix 
B, we get 

6(z > Q(l - z)(M + M N - B x D f N (T- = P- - Hgj^ z, Tl, AQ) 



A D/JV (z, T ± , A ± , Mf, M^)^(z, Tx, -Ax, M 2 , M 2 Di ) 

? A',A / T _ _,p- , (f x -A x /2) 2 + 



< 0(x + 0(M + M N - B x D f N (T~ = £P~ + ^fegS^ , *, Tx, A ± )) 



with 



^ /JV (*, Tx, Ax, Ml M 2 Dl )C D/N (z, Tx, A ± , M^) 

(C.62) 



^(^fx^M^M^) 

(1 - z)(z + 0(1 - 0T 2 - (* + 0(TI + M 2 ) - (1 - z)((fx - Ax/2) 2 + M 2 J, 
C D /iv(^e,Tx,Ax,M 2 J 

20(^ 2 - i 2 )P 2 -Tl- A\/ A - M 2 Di ) - 2zf ± • Ax, 
B x D f N (T-,z,f ± ,A ± )) 



[z(M 2 N - Ai/4) + 2(1 - C)T~P + + £Tx • A ± - t\T ± A AJ 



2M N 

Tx(A) 



+ (Ja'.-aMjv U + 2£ 



Ax(A); ' 
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where P 2 = M 2 N - A 2 /4. 

With help of Eq. (3.26), we finally arrive at 



dy / dz5(x — yz) 



E D (x,£,A 2 ) 



F s D (y,(,Al)-2(F a D (y,CAl) 
G?(A 2 ) 

^ D (y,C,A 2 ,)-2CF f P(y,C,A 2 ,) 



(2^ 

/ dz5(x — yz)- 



(A 2 ) 



(27T) 



where 



9dZn 



E 1 



(*>0(i-*) 



A D /n(z, f , T±, A ± )A D/N (z, 71, -A ± ) 



+ 



1*1 < + 



A D / N (z, f , 71, A±)C D / N (z, f , 71, A^ 



^(T-,z,e,f ± ,A ± ), 



£^(.,£,71,A ± ) 



(*>0(i-*) 



A D/JV (z, £, 71, A ± )A D/JV (z, 71, - A ± ) 



*(M<0(* + 



A D/JV (z, f , 71, A ± )C D/N {z, f , 71, A ± ) 



^(T-,z,e,71,A ± ), 



in which 



H D / N (T-,z,£,f ± ,A ± ) 
zM N 



i-e 



(1 + e)(M + M N ) - C 2 M N (z + 2^f4^) 



A ± (\y 



z{M% - A\l A) + 2(1 - e)T,rP+ + £71 • A 



2M; 



JV 
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with 



Tr . e(z > 0{P - _ J^L, + , N < + ^g±^ ). (c, 7) 

and 

^(7;-, f , Tl, £l) = zMjv (.U • .l/,v - AM* + 2^^) )) • (C- 68 ) 

Note that if we integrate Eq. (C.63) over x with A 2 = A 2 D and fixed £, then we can 
reproduce the diquark current contributions to the form factors, 

JjxH D {x^A 2 ) = J^dzJ^H D / N (z^f ± ,K ± ) = F 1 D (A 2 ), 
J^dxE D (x,£,A 2 ) = £dzf^E D ' N (z,£,f ± ,Z ± ) = F?(A*), 

(C.69) 

where Ff denotes diquark current contributions to the nucleon form factor which are 
given in Appendix E. 
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D Vertex corrections to the photon vertex 

The photon vertex correction, as shown in Fig. 2, consists of the sum of a series of ring 
diagrams. Each diagram on the left side of the diagrams in Fig. 2 is calculated as follows: 

rV + (-2)G a n^(A) abl vT b + (-2) 2 G^n^(A) a6 n^(A) 6c7 ,,r c + • • • , (D.70) 

where 

n^(A) o6 = 6i5 ab J -0- 4 tr D [S(k)^S(k + A)Y\- (D.71) 
Hy(A) ab can be decomposed into the longitudinal and transverse parts: 

ITf (A) o6 = 5 ab [(A 2 g^ - A M A iy )Ily iT (A 2 ) + A^ILv L (A 2 )]. (D.72) 

With the use of Ward identity, then it is clear that the transverse part of Yiy'(A) ab does 
not contribute due to current conservation. Therefore the series can be easily summed: 

(D.73) 



1 + 2G a A 2 riy iT (A 2 ) ' 

where a = 0(i) means the isoscalar (isovector) part and express the coupling constants 
in the vector meson channels. 

Ily i L(A 2 ) can be calculated from 

[(Ay-A"^)n, !T (A 2 )+A"A^ i (A 2 )]=6i J ^L t r D [S{k)Y S{k+A)^}. (D.74) 
Inserting the PV-regularization factor, and performing the /c-integrals, we obtain 



n y , T (A 2 ) = A ( l daail 

71 JO 



a 



A 2 , / A 2 

In 1 + 



Mq - a(l - a) A 2 + A 2 V M 2 - a(l - a)A 2 , 



(D.75) 



From Eqs. (D.73) and (D.75), we can easily see that there is no vertex correction at 
A 2 = 0, i.e., when the photon is on the mass shell. 
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E Nucleon form factors 



Nucleon electromagnetic form factors can be calculated in the same way as the GPDs, with 
the operator 7 + (l ± t z )/2 replaced by 7 M (1 ± t z )/2. We introduce Feynmann parameters 
z, Xi t 2 to combine the denominators of the propagators, and then a Wick rotation is 
performed to obtain an Euclidean integral. The resulting expressions are given by 

pc/a2\ ao r 7 Vr n . i r tdt | + (i-4z 2 )f + M 2 Qi 

F l( A) - -4Q q G s Z N ^ Ci J_ i2 dz-J o ^ [t + Mh _ (1 _ Az2) ^ ? 

if (A 2 ) = -AQ q G s Z N ^ I' dz l - I™ ^ 2 2 2MqMn a , 

21 ' ^ q NZ T J-h 2io 8^[t + M 2 Qi -(l-Az^f}' 

if (A 2 ) = Q q g 2 d Z N ^c t dxj dx 2 - 

; JO J -XI 2 JO 



X 



8tt 2 

(1 - Xl fM 2 N + 2(1 - Xl )M Q M N + M^ + l + (xj - x 2 )f 
[t + (l- Xl )M 2 Dl + x x M 2 Ql - xi(l - Xl )M% - {x 2 - ]s 

toft 

''2 Jo 8^ 

2M jV x 1 [M j y(l- 3 ; 1 ) + MQ] 

[t + (1 - xx)J^ + ^M 2 , - Xl (l - ^M 2 - (x? - x|)^]3 

i^(A 2 ) = q^z^cJ 1 d Xl r dx 2 \r hi! 

; Jo J -xi 2 Jo 



rl rxi 1 rc 

if (A 2 ) = Q q g 2 d Z N Y,Ci dxj dx 2 - 

JO J—xi Z JO 



X 



X 



8tt 2 

2M N (1 - x 1 )(x 1 M N + M Q ) + § 



[t + (l- x x )M 2 Ql + Xl M 2 Dl - Xl (l- x,)M 2 N - (x 2 - xl)f 13 

if (A 2 ) = Q^ZsY^c, f d,- x I" dx 2 \ f 

, Jo J-xi 2 Jo 



x 



-x! '2 Jo 8tt 2 
^M^l-arObiMjv + Mo] 



[t + (1 - Xl )M 2 Qt + x 1 M 2 Di - Xl (l- Xl )M 2 N - (x 2 - x 2 )^} 



where C,P,D mean current, pole and diquark current contributions. 



(E.76) 
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